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Abstract: Recent emergence and growing use of terahertz (THz) radiation 
for medical imaging and public security screening raise questions on 
reasonable levels of exposure and health consequences of this form of 
electromagnetic radiation. In particular, picosecond-duration THz pulses 
have shown promise for novel diagnostic imaging techniques. However, the 
effects of THz pulses on human cells and tissues thus far remain largely 
unknown. We report on the investigation of the biological effects of pulsed 
THz radiation on artificial human skin tissues. We observe that exposure to 
intense THz pulses for ten minutes leads to a significant induction of H2AX 
phosphorylation, indicating that THz pulse irradiation may cause DNA 
damage in exposed skin tissue. At the same time, we find a THz-pulse- 
induced increase in the levels of several proteins responsible for cell-cycle 
regulation and tumor suppression, suggesting that DNA damage repair 
mechanisms are quickly activated. Furthermore, we find that the cellular 
response to pulsed THz radiation is significantly different from that induced 
by exposure to UVA (400 nm). 
©2013 Optical Society of America 

OCIS codes: (170.0170) Medical optics and biotechnology; (170.1420) Biology; (170.6930) 
Tissue; (170.7160) Ultrafast technology. 

References and links 

1. D. M. Mittleman, Sensing with Terahertz Radiation (Springer, 2010). 

2. M. Tonouchi, "Cutting-edge terahertz technology," Nat. Photonics 1(2), 97-105 (2007). 

3. P. U. Jepsen, D. G. Cooke, and M. Koch, "Terahertz spectroscopy and imaging — Modem techniques and 
applications," Laser Photonics Rev. 5(1), 124-166 (201 1). 

4. A. J. Fitzgerald, E. Berry, N. N. Zinovev, G. C. Walker, M. A. Smith, and J. M. Chamberlain, "An introduction 
to medical imaging with coherent terahertz frequency radiation," Phys. Med. Biol. 47(7), R67-R84 (2002). 

5. P. H. Siegel, "Terahertz Technology in Biology and medicine," IEEE Trans. Microw. Theory Tech. 52(10), 
2438-2447 (2004). 

6. G. J. Wilmink and J. E. Grandt, "Current State of Research on Biological Effects of Terahertz Radiation," J. 
Infrared Millim. Terahz Waves 32(10), 1074-1 122 (201 1). 

7. A. G. Markelz, S. Whitmire, J. Hillebrecht, and R. Birge, "THz time domain spectroscopy of biomolecular 
conformational modes," Phys. Med. Biol. 47(21), 3797-3805 (2002). 

8. B. M. Fischer, M. Walther, and P. U. Jepsen, "Far-infrared vibrational modes of DNA components studied by 
teraheilz time-domain spectroscopy," Phys. Med. Biol. 47(21), 3807-3814 (2002). 

9. M. Brucherseifer, M. Nagel, P. Haring Bolivar, H. Kurz, A. Bosserhoff, and R. Biittner, "Label-free probing of 
the binding state of DNA by time-domain terahertz sensing," Appl. Phys. Lett. 77(24), 4049^051 (2000). 

10. S. W. Smye, J. M. Chamberlain, A. J. Fitzgerald, and E. Berry, "The interaction between Terahertz radiation and 
biological tissue," Phys. Med. Biol. 46(9), RlOl-Rl 12 (2001). 

11. M. Nagel, P. Haring Bolivar, M. Brucherseifer, H. Kurz, A. Bosserhoff, and R. Biittner, "Integrated THz 
technology for label-free genetic diagnostics," Appl. Phys. Lett. 80(1), 154-156 (2002). 

12. M. Vasudev, J. Yang, H.-S. Jung, M. A. Stroscio, and M. Dutta, "Integrated nanostracture-semiconductor 
molecular complexes as tools for THz spectral studies of DNA," IEEE Sens. J. 10(3), 524-530 (2010). 



#182102 -$15.00 USD Received 19 Dec 2012; revised 1 Feb 2013; accepted 14 Feb 2013; published 14 Mar 2013 
(C) 2013 OSA 1 April 2013 / Vol. 4, No. 4 / BIOMEDICAL OPTICS EXPRESS 559 



13. S. J. Kim, B. Bom, M. Havenith, and M. Gruebele, "Real-time detection of protein- water dynamics upon protein 
folding by terahertz absorption spectroscopy," Angew. Chem. Int. Ed. Engl. 47(34), 6486-6489 (2008). 

14. A. G. Markelz, "Terahertz Dielectric Sensitivity to Biomolecular Structure and Function," IEEE J. Sel. Top. 
Quantum Electron. 14(1), 180-190 (2008). 

15. P. C. Ashworth, E. Piclcwell-MacPherson, E. Provenzano, S. E. Pinder, A. D. Purushotham, M. Pepper, and V. P. 
Wallace, "Terahertz pulsed spectroscopy of freshly excised human breast cancer," Opt. Express 17( 1 5), 1 2444- 
12454 (2009). 

16. A. J. Fitzgerald, V. P. Wallace, M. Jimenez-Linan, L. Bobrow, R. J. Pye, A. D. Purushotham, and D. D. Amone, 
"Terahertz pulsed imaging of human breast tumors," Radiology 239(2), 533-540 (2006). 

17. R. M. Woodward, V. P. Wallace, D. D. Amone, E. H. Linfield, and M. Pepper, "Terahertz pulsed imaging of 
skin cancer in the time and frequency domain," J. Biol. Phys. 29(2/3), 257-259 (2003). 

18. P. C. Ashworth, P. O'Kelly, A. D. Purashotham, S. E. Pinder, M. Kontos, M. Pepper, and V. P. Wallace, "An 
intra-operative THz probe for use during the surgical removal of breast tumors," in 33rd Intemational Conference 
on Infrared, Millimeter and Terahertz Waves, 2008. IRMMW-THz 2008 (2008), pp. 1-3. doi: 

10.1 109/1C1MW.2008.46658 1 0. 

19. "TeraView receives MHRA approval to conduct in-vivo cancer imaging trials," 
http://www.teraview.com/news/terahertz-news/TeraView-receives-in-vivo-trials-approval.html 

20. J. F. Federici, B. Schulkin, F. Huang, D. Gary, R. Barat, F. Oliveira, and D. Zimdars, "THz imaging and sensing 
for security applications: explosives, weapons and dmgs," Semicond. Sci. Technol. 20(7), S266-S280 (2005). 

21. D. Zimdars and J. S. White, "Terahertz reflection imaging for package and personnel inspection," Proc. SPIE 
5411,78-83 (2004). 

22. B. S. Alexandrov, V. Gelev, A. R. Bishop, A. Usheva, and K. O. Rasmussen, "DNA breathing dynamics in the 
presence of a terahertz field," Phys. Lett. A 374(10), 1214 1217 (2010). 

23. A. Korenstein-llan, A. Barbul, P. Hasin, A. Eliran, A. Gover, and R. Korenstein, "Terahertz radiation increases 
genomic instability in human lymphocytes," Radiat. Res. 170(2), 224-234 (2008). 

24. H. Hintzsche, C. Jastrow, T. Kleine-Ostmann, H. Stopper, E. Schmid, and T. Schrader, "Terahertz radiation 
induces spindle disturbances in human-hamster hybrid cells," Radiat. Res. 175(5), 569-574 (201 1). 

25. G. J. Wilmink, B. L. Ibey, C. L. Roth, R. L. Vincelette, B. D. Rivest, C. B. Hom, J. Bemhard, D. Roberson, and 
W. P. Roach, "Determination of death thresholds and identification of terahertz (THz) - specific gene expression 
signatures," Proc. SPIE 7562, 75620K, 75620K-8 (2010). 

26. G. J. Wilmink, B. D. Rivest, B. L. Ibey, C. L. Roth, J. Bemhard, and W. P. Roach, "Quantitative investigation of 
bioeffects associated with terahertz radiation," Proc. SPIE 7562, 75620L, 75620L-10 (2010). 

27. J. Bock, Y. Fukuyo, S. Kang, M. L. Phipps, L. B. Alexandrov, K. 0. Rasmussen, A. R. Bishop, E. D. Rosen, J. S. 
Martinez, H.-T. Chen, G. Rodriguez, B. S. Alexandrov, and A. Usheva, "Mammalian stem cells reprogramming 
in response to terahertz radiation," PLoS ONE 5(12), el5806 (2010). 

28. B. S. Alexandrov, K. 0. Rasmussen, A. R. Bishop, A. Usheva, L. B. Alexandrov, S. Chong, Y. Dagon, L. G. 
Booshehri, C. H. Mielke, M. L. Phipps, J. S. Martinez, H.-T. Chen, and G. Rodriguez, "Non-thermal effects of 
terahertz radiation on gene expression in mouse stem cells," Biomed. Opt. Express 2(9), 2679-2689 (201 1). 

29. E. Berry, G. C. Walker, A. J. Fitzgerald, N. N. Zinov'ev, M. Chamberlain, S. W. Smye, R. E. Miles, and M. A. 
Smith, "Do in vivo terahertz imaging systems comply with safety guidelines?" J. Laser Appl. 15(3), 192-198 
(2003). 

30. T. T. L. Kristensen, W. Withayachumnankul, P. U. Jepsen, and D. Abbott, "Modeling terahertz heating effects on 
water," Opt. Express 18(5), 4727-4739 (2010). 

31. S. M. Chitanvis, "Can low-power electromagnetic radiation dismpt hydrogen bonds in dsDNA?" J. Polym. Sci., 
B, Polym. Phys. 44(18), 2740-2747 (2006). 

32. W. Zhuang, Y. Feng, and E. W. Prohofsky, "Self-consistent calculation of localized DNA vibrational properties 
at a double-helix-single-sttand junction with anharmonic potential," Phys. Rev. A 41(12), 7033-7042 (1990). 

33. E. P. Rogakou, D. R. Pilch, A. H. Orr, V. S. Ivanova, and W. M. Bonner, "DNA double-stranded breaks induce 
histone H2AX phosphorylation on serine 139," J. Biol. Chem. 273(10), 5858-5868 (1998). 

34. O. A. Sedelnikova, D. R. Pilch, C. Redon, and W. M. Boimer, "Histone H2AX in DNA damage and repair," 
Cancer Biol. Ther. 2(3), 233-235 (2003). 

35. F. Blanchard, G. Sharma, L. Razzari, X. Ropagnol, H.-C. Bandulet, F. Vidal, R. Morandotti, J.-C. Kieffer, T. 
Ozaki, H. Tiedje, H. Haugen, M. Reid, and F. A. Hegmann, "Generation of intense terahertz radiation via optical 
methods," IEEE J. Sel. Top. Quantum Electron. 17(1), 5-16 (201 1). 

36. M. C. Hoffmann and J. A. Fiilop, "Intense ulttashort terahertz pulses: generation and applications," J. Phys. D 
Appl. Phys. 44(8), 083001 (201 1). 

37. E. Boelsma, S. Gibbs, C. Faller, and M. Ponec, "Characterization and comparison of reconstmcted skin models: 
morphological and immunohistochemical evaluation," Acta Derm. Venereol. 80(2), 82-88 (2000). 

38. O. A. Sedelnikova, A. Nakamura, O. Kovalchuk, 1. Koturbash, S. A. Mitchell, S. A. Marino, D. J. Brenner, and 
W. M. Bonner, "DNA double-strand breaks form in bystander cells after microbeam irradiation of three- 
dimensional human tissue models," Cancer Res. 67(9), 4295^302 (2007). 

39. N. S. Agar, G. M. Halliday, R. S. Bametson, H. N. Ananthaswaray, M. Wheeler, and A. M. Jones, "The basal 
layer in human squamous tumors harbors more UVA than UVB fingerprint mutations: a role for UVA in human 
skin carcinogenesis," Proc. Natl. Acad. Sci. U.S.A. 101(14), 4954-4959 (2004). 



#182102 - $15.00 USD Received 19 Dec 2012; revised 1 Feb 2013; accepted 14 Feb 2013; published 14 Mar 2013 
(C) 2013 OSA 1 April 2013 / Vol. 4, No. 4 / BIOMEDICAL OPTICS EXPRESS 560 



40. S. Q. Wang, R. Setlow, M. Berwick, D. Polsky, A. A. Marghoob, A. W. Kopf, and R. S. Bart, "Ultraviolet A and 
melanoma: a review," J. Am. Acad. Dermatol. 44(5), 837-846 (2001). 

41. J. Cadet, E. Sage, and T. Douki, "Ultraviolet radiation-mediated damage to cellular DNA," Mutat. Res. 571(1-2), 
3-17(2005). 

42. R. P. Rastogi, A. Richa, A. Kumar, M. B. Tyagi, and R. P. Sinha, "Molecular mechanisms of ultraviolet 
radiation-induced DNA damage and repair," J. Nucleic Acids 2010, 592980 (2010). 

43. O. Kovalchuk, F. J. Zemp, J. N. Filkowski, A. M. Altamirano, J. S. Dickey, G. Jenkins-Baker, S. A. Marino, D. J. 
Brenner, W. M. Bonner, and O. A. Sedelnikova, "microRNAome changes in bystander three-dimensional human 

tissue models suggest priming of apoptotic pathways," Carcinogenesis 31(10), 1882-1888 (2010). 

44. J. H. J. Hoeijmakers, "Genome maintenance mechanisms for preventing cancer," Nature 411(6835), 366-374 
(2001). 

45. J. S. Dickey, F. J. Zemp, A. Altamirano, O. A. Sedelnikova, W. M. Bonner, and O. Kovalchuk, "H2AX 
phosphorylation in response to DNA double-strand break formation during bystander signalling: effect of 
microRNA knockdown," Radiat. Prot. Dosimetry 143(2-4), 264-269 (2011). 

46. W. M. Bonner, C. E. Redon, J. S. Dickey, A. J. Nakamura, O. A. Sedelnikova, S. Solier, and Y. Pommier, 
"yH2AX and cancer," Nat. Rev. Cancer 8(12), 957-967 (2008). 

47. H. T. Wang, B. Choi, and M. S. Tang, "Melanocytes are deficient in repair of oxidative DNA damage and UV- 
induced photoproducts," Proc. Natl. Acad. Sci. U.S.A. 107(27), 12180-12185 (2010). 

48. J. Fillingham, M.-C. Keogh, and N. J. Krogan, "yH2AX and its role in DNA double-strand break repair," 
Biochem. Cell Biol. 84(4), 568-577 (2006). 

49. P. McGlynn and R. G. Lloyd, "Recombinational repair and restart of damaged replication forks," Nat. Rev. Mol. 
Cell Biol. 3(1 1), 859-870 (2002). 

50. K. K. Khanna and S. P. Jackson, "DNA double-strand breaks: signaling, repair and the cancer connection," Nat. 
Genet. 27(3), 247-254 (2001). 

51. D. Menendez, A. Inga, and M. A. Resnick, "The expanding universe of p53 targets," Nat. Rev. Cancer 9(10), 
724-737 (2009). 

52. E. Bolderson, D. J. Richard, B. B. Zhou, and K. K. Khanna, "Recent advances in cancer therapy targeting 
proteins involved in DNA double-strand break repair," Clin. Cancer Res. 15(20), 63 14-6320 (2009). 

53. S. Seite, D. Moyal, M. P. Verdier, C. Hourseau, and A. Fourtanier, "Accumulated p53 protein and UVA 
protection level of sunscreens," Photodermatol. Photoimmunol. Photomed. 16(1), 3-9 (2000). 

54. A. Krones-Herzig, S. Mittal, K. Yule, H. Liang, C. English, R. Urcis, T. Soni, E. D. Adamson, and D. Mercola, 
"Early growth response 1 acts as a tumor suppressor in vivo and in vitro via regulation of p53," Cancer Res. 
65(12), 5133-5143 (2005). 

55. Y. Zwang, A. Sas-Chen, Y. Drier, T. Shay, R. Avraham, M. Lauriola, E. Shema, E. Lidor-Nili, J. Jacob-Hirsch, 
N. Amariglio, Y. Lu, G. B. Mills, G. Rechavi, M. Oren, E. Domany, and Y. Yarden, "Two phases of mitogenic 
signaling unveil roles for p53 and EGRl in elimination of inconsistent growth signals," Mol. Cell 42(4), 524-535 
(2011). 

56. B. L. Mahaney, K. Meek, and S. P. Lees-Miller, "Repair of ionizing radiation-induced DNA double-strand 
breaks by non-homologous end-joining," Biochem. J. 417(3), 639-650 (2009). 

57. P. Reynolds, J. A. Anderson, J. V. Harper, M. A. Hill, S. W. Botchway, A. W. Parker, and P. O'Neill, "The 
dynamics of Ku70/80 and DNA-PKcs at DSBs induced by ionizing radiation is dependent on the complexity of 
damage," Nucleic Acids Res. 40(21), 10821-10831 (2012). 

58. H. Chen, Y. Bao, L. Yu, R. Jia, W. Cheng, and C. Shao, "Comparison of cellular damage response to low-dose- 
rate 1251 seed irradiation and high-dose-rate gamma irradiation in human lung cancer cells," Brachytherapy 
11(2), 149-156 (2012). 

59. C. J. Sherr and J. M. Roberts, "Inhibitors of mammalian Gl cyclin-dependent kinases," Genes Dev. 9(10), 1 149- 
1163 (1995). 

60. T. Sperka, J. Wang, and K. L. Rudolph, "DNA damage checkpoints in stem cells, ageing and cancer," Nat. Rev. 
MoL Cell Biol. 13(9), 579-590 (2012). 



1. Introduction 

Recent technological breakthroughs have led to numerous applications of THz radiation [1- 
6] . In the medical field, sensitivity of THz radiation to the molecular properties of biological 
tissue, from cell hydration to conformation states of important biomolecules [7-14], has 
spurred development of new THz-based approaches to cancer diagnosis. AppHcation of THz 
pulsed imaging to detection of breast cancer [15], including in situ breast cancer [16], as well 
as melanoma [4] and basal cell skin carcinoma [17] has been demonstrated. The first in vivo 
trials of THz imaging as an intra-operative tool during cancer surgery are currently underway 
[18,19]. In security applications, THz radiation is now extensively applied for identification of 
concealed explosives and weapons [20,21]. While these new THz radiation-based applications 
have generated excitement in the research community, concerns have been raised regarding 
the possible health risks associated with THz exposure [22-28]. 
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The current understanding of biological effects of THz radiation, i.e., its potential to 
induce DNA damage and impact cell activity, is still limited. Until recently, it was assumed 
that exposure risks are thermal in nature due to absorption of THz radiation by water in 
biological tissue [4,6,29,30]. However, several theoretical studies have suggested that while 
the energy of THz photons is too low to break chemical bonds, resonance-type linear and 
nonlinear interactions of THz electromagnetic fields with DNA may, under certain conditions, 
significantly alter DNA dynamics and even induce localized openings (bubbles) in the DNA 
strands [22,31,32]. This is of particular concern for applications of intense THz pulses of 
picosecond duration. A few experimental studies have explored the cellular and molecular 
responses to continuous-wave, monochromatic THz radiation and found evidence of THz- 
induced genotoxicity in human lymphocytes [23], spindle disturbances in human-hamster 
hybrid cells [24], as well as changes in gene expression and activation of apoptotic and 
necrotic processes in human dermal fibroblasts and Jurkat cells [25,26]. A recent study 
exposing mouse stem cells to broadband THz pulses reported the first experimental 
confirmation of THz-pulse-induced gene expression changes in mammalian cells, possibly 
correlated to THz-induced breathing vibrational modes in the corresponding promoter DNA 
[27,28]. However, the effects of intense broadband (picosecond) THz pulses on human cells 
and tissues are not known. In particular, since the penetration depth of THz radiation into the 
human body is limited to a fraction of a millimeter [5], knowledge of the biological effects of 
intense THz pulses used in novel cancer detection modalities on human skin is important. 

In this work, we show evidence strongly indicative of double strand breaks (DSBs) in 
DNA induced by intense, picosecond THz pulses in exposed artificial human skin tissue 
models. We use the presence of phosphorylated H2AX (yH2AX), which is one of the earliest 
and most characterized cellular responses to DSBs [33,34], as a surrogate marker for DNA 
damage. At the same time, we observe THz-pulse-induced increases in the levels of multiple 
tumor suppressor and cell-cycle regulatory proteins that facilitate DNA repair. This may 
suggest that DNA damage in human skin arising from broadband THz pulse exposure could 
be quickly and efficiently repaired, therefore minimizing the risk of point mutation, prelude to 
carcinogenesis. We note that the energy of the intense THz pulses used in this study of 0.1 - 

1.0 (xJ is many orders of magnitude higher than pulse energies of approximately 10"** |aJ used 
in current medical imaging applications that use traditional photoconductive THz pulse 
emitters [30]. However, the observed capability of intense THz pulses to affect DNA and 
cellular functions warrants exploration into potential therapeutic apphcations of intense THz 
pulses. 

2. Materials and methods 

2.1 THz pulse source, tissue models, and exposure conditions 

In our experiments, we exposed artificial human 3D skin tissues to broadband THz pulses 
with -1.7 ps duration, 1 kHz repetition rate, and pulse energy variable up to 1 \iS. The THz 
pulses were generated by optical rectification of tilted-pulse-front 800 nm pulses from an 
amplified Ti:sapphire laser source in LiNbOs [35,36], as shown in Fig. 1(a). A pyroelectric 
detector (Spectrum Detector) was used to measure the THz pulse energies. Free-space electro- 
optic sampling was used to record the THz pulse waveform, as shown in Fig. 1(b). The 
corresponding amplitude spectrum of the THz pulse, which is shown in Fig. 1(c), is peaked at 
0.5 THz and has a bandwidth of 0.1 - 2 THz. The full thickness human skin tissue model 
(EpiDermFT, MatTek) reconstructs normal tissue structure and consists of normal, human- 
derived epidermal keratinocytes and dermal fibroblasts that form a multilayered, highly 
differentiated model of human dermis and epidermis, as shown in Fig. 1(d). It is mitotically 
and metabolically active, and preserves the arrangement and communication of cells in skin 
tissue in vivo [37,38], and thus provides an appropriate platform for assessing the effects of 
exposure to high intensity THz pulses on human skin. The tissues were cultured according to 
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the manufacturer's protocol, using a conventional air-liquid interface tissue culture technique. 
The top surface of the tissue was exposed to air, and the bottom surface of the tissue was fed 
by a culture medium through pores in a membrane. Tissue samples in single well plate inserts 
were removed from the medium and placed at the focus of the 1 kHz pulsed THz beam with 
pulse energies of either 1 \xJ (referred to here as "high" THz energy) or 0.1 \xJ ("low" THz 
energy) for 10 minutes. The THz radiation was incident on the top surface of the tissue 
sample. The THz spot size at the focus was 1.5 mm diameter, as imaged by a pyroelectric 
infrared camera (Electrophysics PV320) and shown in Fig. 1(e). This produced peak incident 
THz electric fields at the sample of approximately 220 kV/cm and 70 kV/cm at high and low 
THz pulse energies, respectively. Since the penetration depth of THz radiation in human skin 
is less than 150 |am [5], most of the THz pulse energy (> 93%) is absorbed within the 400 nm 
thick tissue samples. 




Time (ps) Frequency (THz) 





_ Stratum 
corneum 

Epidermis (e) 



■ Dermis 

100 pm ■ '• 



Fig. 1. Exposure of artificial human skin tissues equivalents to intense THz pulses, (a) 
Schematic of the tilted-pulse-front THz pulse source and EpiDermFT tissue in a single-well 
inseil placed at the focus of THz beam, (b) Waveform of a 1.0 pJ THz pulse and (c) 
coiTcsponding amplitude spectrum, (d) Histology of the EpiDermFT tissue sample (400X 
image by MatTek Corporation, used with permission). The stained paraffin section reveals 
epidermis containing basal, spinous, granular keratinocytes and stratum corneum. The dermis 
contains numerous viable fibroblasts, (e) Image of the THz beam spot size at the sample 
location. The 1/e^ diameter of the THz beam is 1.5 mm. The dashed square shows the size of 
the area (2x2 mm^) that was cut out to separate exposed from unexposed tissue. 

Exposure to UVA, a known genotoxic agent, was performed using 400 nm, 50 fs, 0.080 i^J 
pulses with 1 kHz repetition rate for 2 minutes and was used as a positive control for 
induction of DNA damage [39-42]. The UVA beam spot size at the tissue surface was 2.7 
mm diameter. Unexposed tissue samples served as controls (CT). Four tissue samples were 
used for each of the four experimental conditions (CT, Low, High, and UVA). The tissues 
were then stored in a multiwell dish filled with fresh medium and incubated at 37 °C for 30 
minutes, followed by a snap freeze on dry ice. Control tissues underwent the same procedures 
other than being irradiated. THz exposures were carried out at biologically-low ambient 
temperature (21 °C). The time-averaged THz power density was 57 mW/cm^ at the highest 
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THz pulse energy. Applying the theoretical formalism developed in [30], we estimate that 
temperature increase due to THz exposure was less than 0.7 °C. 

2.2 Protein extraction and western immunoblotting 

The central square region (2 mm x 2 mm, Fig. 1(e)) containing the exposed portion of the 
frozen tissues was cut out and used for western immunoblotting analysis for determining the 
levels of various cell-cycle regulatory and tumor suppressant proteins, as described in [43]. 
Briefly, tissue samples were sonicated in 0.4 ml of ice-chilled 1% sodium dodecyl sulphate 
(SDS) and boiled for 10 min. Small aliquots (10 |il) of homogenate were reserved for protein 
determination using assay reagents from BioRad (Hercules, CA, USA). Since the amount of 
starting tissue was very small, upon quantification protein extracts of three independent 
biological replicates were pooled together with equal representation of each sample. Pooled 
protein extracts were separated by SDS-polyacrylamide electrophoresis (PAGE) in slab gels 
of 12% polyacrylamide, made in duplicates, and transferred to PVDF membranes (Amersham, 
Bale d'Urfe, QC, Canada). Membranes were incubated with antibodies against yH2AX, p21, 
pl6, pl5, p27, EGRl (1:1000, Cell Signaling, Danvers, MA, USA), p53 (1:500, Santa Cruz 
Biotechnology, Santa Cruz, CA), p-p53 (Serl5) (1:500, Cell Signaling) KU70 (1:1000, Santa 
Cruz Biotechnology), actin (1:2000, Abeam, Cambridge, MA). Antibody binding was 
revealed by incubation with horseradish peroxidase-conjugated secondary antibodies 
(Amersham) and ECL Plus immunoblotting detection (Amersham). Chemiluminescence was 
detected by a FluorChem HD2 System (Cell Biosciences/ProteinSimple, Santa Clara, CA, 
USA) and Biomax MR films (Eastman Kodak, New Haven, CT, USA). Unaltered PVDF 
membranes were stained with Coomassie Blue (BioRad). 

2.3 y-H2AX immunochemistry 

For each distinctive exposure type condition (high, low, UVA, CT), one frozen tissue was 
sectioned into 5-10-nm-thick slices for immunochemistry analysis, as described in [38]. 
Briefly, for the y-H2AX assay, fi-ozen sections were dried, fixed in 2% paraformaldehyde, 
permeabilized with 1% Triton X-100, and processed for immunostaining with an anti-yH2AX 
rabbit antibody. The nuclei were counterstained with 4,6-diamidrno-2-phenylindole (DAPI). 

3. Results and discussion 

Earlier theoretical studies have suggested that exposure to high intensity THz radiation with 
frequencies resonant with the local breathing dynamics of double-stranded DNA may lead to 
the formation of localized openings in the double strands [22]. Picosecond duration, intense 
THz pulses have instantaneous THz fields of hundreds of kV/cm, which might be sufficient to 
not only temporarily disrupt DNA dynamics but also cause DSBs to occur. In DSB, two 
complementary strands of the DNA double helix are damaged simultaneously. It is potentially 
the most dangerous type of DNA damage, which, if unrepaired, can lead to cell death or 
cancer [34]. DSBs can occur in cells due to many different endogenous factors, such as 
oxidative damage by reactive oxygen species, as well as exogenous factors. Of the known 
exogenous factors, ionizing radiation, such as X-rays [44] or UV radiation [39-41] are the 
most prominent causes of DSBs. Regardless of the exact mechanism of DSB formation, they 
can be experimentally detected by assaying for the presence of the phosphorylated histone 
H2AX (yH2AX). Hundreds of H2AX molecules in the chromatin surrounding the break site 
become phosphorylated at the Serine 139 position within minutes of DSB formation [34,45]. 
The resulting yH2AX foci attract repair proteins to the DSB site, thus playing a critical role in 
DSB recognition and initiation of repair processes [34]. Analysis of yH2AX expression is 
widely used as a surrogate marker for DNA damage and a reliable DSB detection tool, since it 
is orders of magnitude more sensitive compared to other DSB visualization techniques 
[34,46]. 
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To assess the occurrence of DSB in THz-exposed and UV-exposed tissues compared to 
unexposed controls, we assayed for the presence of yH2AX using western blot analysis. We 
observe that exposure to intense THz pulses for just ten minutes leads to a significant 
induction of H2AX phosphorylation compared to control samples, as seen in Western blot 
images for both high and low THz pulse energy exposures in Fig. 2. This is an indication that 
intense THz pulses cause DNA DSBs to occur in human skin cells. UVA (400 nm), which 
creates reactive oxygen species and causes oxidative DNA damage [40,41,47], also leads to 
significant increases in yH2AX levels in exposed tissues (Fig. 2). 
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Fig. 2. Induction of yH2AX in artificial human skin tissues equivalents following 10 minute 
exposure to either high (1.0 or low (0.1 energy THz pulses, or 2 minute exposure to 
UVA (400 nm) pulses (0.080 as compared to control (CT) samples. The pulse repetition 
rate was 1 kHz. Actin was used as a loading control. Each experiment included pooled lysates 
from three tissues for each exposure condition with equal representation for each tissue. 

The presence and localization of DSB sites in exposed tissue samples can be also studied 
by visualizing y-H2AX-expressing cells by in situ immunofluorescence with yH2AX 
antibodies. Figure 3 shows representative images of CT, high, low, and UVA sample cross- 
sections. In these images, yH2AX foci appear as bright green, and DNA as bright blue spots. 
Co-localization of green (yH2AX) and blue (DNA) staining indicates sites within the cells 
where DNA DSBs likely occurred. In UVA-exposed tissue, we observed DNA damage in the 
uppermost epidermis, as expected [39-41]. In tissue exposed to high energy THz pulses, 
multiple yH2AX expressing cells appear in the top layer of epidermis. In low THz pulse 
energy-exposed tissue, yH2AX expressing cells are seen in the top epidermal layer, as well as 
in the dermis. As only one tissue sample per each exposure condition was dedicated to 
immunochemical analysis, we are unable to carry out rehable statistical analysis of the 
incidence of yH2AX-expressng cells within the tissues. While more systematic 
immunochemical analysis of THz-exposed skin tissues is necessary to get a quantitative 
measure of yH2AX expressing cells within the tissues, their depth distribution and exposure 
parameter dependence, this data set serves as an independent verification of THz-pulse 
induced DNA damage in artificial human skin tissues. 
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Fig. 3. Confocal fluorescence images of tissue sample cross-sections showing the incidence of 
7H2AX foci (localized green spots). Blue spots indicate cell nuclei. yH2AX-harboring cells 
(markers for double strand breaks in DNA) inside cell nuclei appear as blue-green spots. Panels 
from left to right show representative images for control (CT), low THz pulse energy (0.1 |.iJ), 
high THz pulse energy (1.0 ^iJ), and UVA-exposed tissues. Scale bar is 100 nm. 
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Having seen the THz-induced increase of yH2AX levels, which are indicative of DNA 
damage, we next proceeded to analyze possible impacts on cellular function by assaying the 
levels of proteins involved in cell-cycle control and DNA repair. Numerous mechanisms exist 
by which cells can repair DSBs. Once DNA damage occurs, checkpoints that pause the cell 
cycle and allow time for DNA repair are activated [46,48]. Two pathways by which DSBs can 
be repaired are homologous recombination and nonhomologous end joining, as reviewed in 
[44,49,50]. The latter mechanism is faster and dominates DSB repair, but is susceptible to 
error as it links broken DNA ends together without reference to accurate base pairing [44]. In 
order to infer how efficient repair of DNA damage caused by intense THz pulse exposure may 
be, we analyzed how THz exposure affects the expression level of various proteins that take 
part in cell cycle control and DSB repair processes by western immunoblotting (Fig. 4). 
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Fig. 4. Expression of tumor suppressor and cell cycle regulatory proteins in artificial human 3D 
skin tissues equivalents following 10 minute exposure to either high (1.0 |xJ) or low (0.1 nJ) 
energy THz pulses, or two minute exposure to 0.080 (jJ UVA (400 nm) pulses, as compared to 
unexposed control (CT) samples. Actin was used as a loading control. 



One of our most striking findings is an increase in the levels of protein p53, which is 
significantly more pronounced in THz-exposed tissues compared to UVA-exposed tissues 
(Fig. 4). An important tumor suppressor and cell-cycle regulator, p53 has been termed 'the 
guardian of the genome' [51,52]. It regulates the cell growth arrest and apoptosis, and also 
takes part in DSB repair [51,52]. Previous studies have reported that repetitive exposures to 
continuous-wave UVA induced a strong increase in the levels of p53-positive cells in the 
human epidermis after a total UVA dose of over 12.5 J/cm^ [53], 100 fold larger than our two 
minute exposure to pulsed UVA hght, which we observed was insufficient to cause significant 
upregulation of p53. On the other hand, the level of the phosphorylated variant of the p53 
protein, p-p53, is barely affected in THz-exposed tissues compared to controls, as shown in 
Fig. 4, and significantly increases with UV exposure. Here, we measured the levels of p53 
phosphorylation of Serine 15 (p-p53 Serl5). This modification is usually catalysed by 
damage-activated ATM kinase and further transactivates the kinase inhibitor p21, leading to 
cell-cycle arrest at the Gl/S checkpoint. Moreover, p53 may influence the G2/M checkpoint 
via transactivation of p21. In our experiment, the tissues were frozen 30 minutes after THz 
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exposure, which might not be sufficient time for p53 phosphorylation at Serl5 by ATM. 
Alternatively, THz-induced DNA damage may not activate ATM, or may cause its delayed 
activation. Another potent in vivo tumor suppressor and p53 target, early growth response 1 
(EGRl) [54], is also slightly upregulated by THz exposure, albeit weaker than by UVA. 
Recent studies have shown that p53 and EGRl collaborate in controlling cellular growth and 
cell cycle, thus allowing time for DNA repair [55]. In mammals, most DSBs are repaired by 
nonhomologous end joining [44,50,56]. Protein KU70 is essential for this process and is 
induced in mammalian cells exposed to X-rays [57], gamma-rays [58], as well as UV 
radiation [42]. We note a significant upregulation of KU70 in THz and UVA-exposed tissues, 
as shown in Fig. 4, indicative of activation of nonhomologous end joining. 

We have also measured the levels of four key cell cycle regulator proteins, p21, pl6, pl5 
and p27, in THz and UVA exposed tissues. These all belong to the cyclin-dependent kinase 
inhibitor protein family, and act as tumor suppressors [59]. Their upregulation provides cells 
crucial time for repairing DNA. We note an increase in the levels of p21, pl6, and p27 in 
THz-exposed tissues but not in UVA-exposed tissues. The levels of another cell cycle 
regulatory protein, pi 5, were unaltered by either THz pulses or UVA. The observed concerted 
upregulation of p53 (discussed earlier), pl6 and p21 is a hallmark of DNA damage response 
in cells [60]. Thus, the observed simultaneous increase in expression of multiple known tumor 
suppressor and damage repair proteins (p53, p21, pl6, p27 and KU70) by intense THz pulse 
exposure suggests that DNA damage was induced, and that damage repair mechanisms were 
promptly activated. However, while further study is required to elucidate the specific roles of 
the aforementioned proteins in cellular response to THz pulse exposure, the observed ability 
of intense THz pulses to induce DNA damage and affect cellular repair functions naturally 
leads to the question; Can THz pulses be used to treat of diseases such as cancer? Given that 
the penetration depth of THz waves into the human body is sub-mm, exploring the possibility 
of the application of THz pulses to treat skin cancers holds the most promise for having 
significant impact on patient care. 

Furthermore, it is known that UVA radiation damages DNA via photosensitized reactions 
that form reactive oxygen species [40-42]. While this damage is repaired by base excision, it 
often yields abasic sites and single-strand breaks, which may later be converted to DSBs 
during DNA replication [42]. THz photons have energies three orders of magnitude smaller 
than UVA photons, which is insufficient to directly break covalent bonds and form free 
radicals. However, recent studies have shovm that DNA, proteins, as well as other important 
biomolecules have collective structural vibrational modes at THz frequencies [7,8,10-14]. 
Resonant coupling of intense THz pulses to those intermolecular vibrational modes may 
therefore induce new conformation states of various cellular biomolecules, or destabihze the 
double-stranded DNA structure via nonlinear resonant effects as suggested in [22], thereby 
affecting cellular functions. More theoretical and experimental studies are clearly needed to 
elucidate the mechanisms by which THz pulses interact with living cells and tissues. 

4. Summary 

In summary, we have observed that exposure to intense THz pulses induces phosphorylation 
of H2AX, indicative of the formation of DNA double strand breaks, and at the same time 
profoundly activates DNA damage response in artificial human skin tissues. Simultaneous 
upregulation of multiple important tumor suppressor proteins in the exposed skin strongly 
suggests that THz-pulse-induced DNA damage may be effectively repaired. Future studies are 
needed to understand the precise molecular mechanisms for the induction of DNA damage 
and changes in cellular signaling by intense picosecond THz pulses, as well as the cellular 
repercussions of such changes. This study may therefore serve as an important roadmap for 
future analysis of the mechanisms and cellular effects of pulsed THz exposure in human 
tissue. 
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